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We show that non-conventional statistical effects (due to the presence of long range 
forces, memory effects, correlations and fluctuations) can be very relevant in the 
interpretation of the experimental observables in relativistic heavy-ions collisions. 
Transverse mass spectrum, transverse momentum fluctuations and rapidity spectra 
are analysed in the framework of the non-extensive statistical mechanics. 



1 Introduction and motivations 

Most of the theoretical analyses related to observables in relativistic heavy- 
ion collisions involve (implicitly or explicitly) the validity of the standard 
Boltzmann-Gibbs statistical mechanics. In particular, if the thermal equilib- 
rium is achieved, the Maxwell-Boltzmann (MB) distribution (Fermi-Dirac or 
Bose-Einstein distribution if quantum statistical effects are not negligible) is 
assumed to hold. When the system approaches equilibrium, the phase-space 
distribution should be derived as a stationary state of the dynamical kinetic 
evolution equation. It is well known that in the absence of non-Markovian 
memory effects, long-range interactions and local correlations, the MB distri- 
bution is obtained as a steady state solution of the kinetic Boltzmann equa- 
tion. However, it is a rather common opinion that, because of the extreme 
conditions of density and temperature in ultrarelativistic heavy ion collisions, 
memory effects and long-range color interactions give rise to the presence of 
non-Markovian processes in the kinetic equation affecting the thcrmalization 
ggggss toward equilibrium as well as the standard equilibrium distribution 

The aim of the present contribution is to explore, from a phcnomenological 
point of view, the relevance of the above mentioned statistical effects that can 
influence the dynamical evolution of the generated fireball toward the freeze- 
out stage and, as a consequence, the physical observables. 
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2 Generalized non extensive statistics 



A quite interesting generalization of the conventional Boltzmann-Gibbs statis- 
tics has been recently proposed by Tsallis □ and proves to be able to overcome 
the shortcomings of the conventional statistical mechanics in many physical 
problems, where the presence of long-range interactions, long-range micro- 
scopic memory, or fractal space-time constraints hinders the usual statistical 
assumptions. 



The Tsallis generalize! 
cralization of the entropy I 



thermostatistics is based upon the following gen- 



w 



s 9 = -f I XX 1 -P? )' W 

where pi is the probability of a given microstate among W different ones and 
q is a fixed real parameter. The new entropy has the usual properties of 
positivity, equiprobability, concavity and irreversibility, preserves the whole 
mathematical structure of thermodynamics and reduces to the conventional 
Boltzmann-Gibbs entropy S = — ^jPilogpi in the limit q — ► 1. 

The single particle distribution function is obtained through the usual 
procedure of maximizing the Tsallis entropy under the constraints of keeping 
constant the average internal energy and the average number of particles. For 
a dilute gas of particles and/or for g w 1 values, the average occupational 
number can be written in a simple analytical form 

{ni)q = [l+tg-W-^-Dll ' (2) 

where the + sign is for fermions, the — for bosons and /3 = 1/T. In the limit 
q — > 1 (extensive statistics), one recovers the conventional Fermi-Dirac and 
Bose-Einstein distribution. Under the same conditions, but in the classical 
limit, one has the following generalized Maxwell-Boltzmann distribution 



^^[l + iq-lWEi- »)] 1/{1 - q) • (3) 
When the entropic q parameter is smaller than 1, the distributions '(Q) and 
(||) have a natural high energy cut-off: Ei < l/[/3(l — q)] + fi, which implies 
that the energy tail is depleted; when q is greater than 1, the cut-off is absent 
and the energy tail of the particle distribution (for fermions and bosons) is 
enhanced. Hence the nonextensive statistics entails a sensible difference of the 
particle distribution shape in the high energy region with respect to the stan- 
dard statistics. This property plays an important role in the interpretation of 
the physical observables, as it will be shown in the following. 
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3 Transverse mass spectrum and momentum fluctuations 



Let us consider the transverse momentum distribution of particles produced, 
e.g., in relativistic heavy ion collisions: it depends on the phase-space distri- 
bution and usually an exponential shape is employed to fit the experimental 
data. This shape is obtained by assuming a purely thermal source with a MB 
distribution. High energy deviations from the exponential shape are taken 
into account by introducing a dynamical effect due to collective transverse 
flow, also called blue-shift. 

Let us consider a different point of view and argue that if long tail time 
memory and long-range interactions are present, the MB distribution must 
be replaced by the generalized distribution (||) . Limiting ourselves to consider 
here only small deviations from standard statistics (q — 1 w 0); then at first 
order in (g — 1) the transverse mass spectrum can be written as 

^ = C -H Z) + ^ P *(*) + *,(*)]} , (4) 

where Ki are the modified Bessel function at the i-order. 

The above equation is able to reproduce very well the transverse mo- 
mentum distribution of hadcpns produced in S+S collisions (NA35 data □) 
providing we take q — I .038 0. Furthermore, it is easy to see that at first or- 
der in (q— I) from Eq.(||), the generalized slope parameter takes the following 
form: 

T q = T + (q-l)m ± . (5) 

Hence nonextensive statistics predicts, in a purely thermal source, a general- 
ized g-blue shift factor at high m± ; moreover this shift factor is not constant 
but increases (if q > I) with m± = m 2 + p\, where m is the mass of the 
detected particle. Such a behavior has been observed in the experimental 
NA44 results!. 

Another observable very sensitive to non-conventional statistical effects 
are the transverse momentum fluctuations. In fact, in the framework of non- 
extensive statistics, the particle fluctuation (An 2 } q = (n 2 ) q — (n) 2 is deformed, 
with respect to the standard expression, as follows 

/A 2\ _ 1 d ( n )q (">9 / -i _j_ / \ \ ra\ 

^ ^ = = l + ( q - m E-,l) (1 ' (6) 



where E is the relativistic energy E = y/ m 2 + p 2 . 

The fluctuations of an ideal gas of fermions (bosons), expressed by Eq.(^), 
are still suppressed (enhanced) by the factor 1 =F (n) q (as in the standard case) 
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but this effect is modulated by the factor [1 + (q — l)f3(E — A*)] -1 - Therefore 
the fluctuations turn out to be increased for q < 1 and are decreased for q > 1 . 
Very good agreement with the experimental NA49 analysis El is obtained by 
taking q = 1.038 cl, notably the same value used in transverse momentum 
spectra. 

4 Anomalous diffusion in rapidity spectra 

An important observable in relativistic heavy-ion collisions is the rapidity 
distribution of the detected particles. In particular, there is experimental 
and theoretical evidence that the broad rapidity distribution of net proton 
yield (p — p) in central heavy-ion collisions at SPS energies could be a signal 
of non-equilibrium properties of the system. We want to show now that 
the broad rapidity shape can be well reproduced in the framework of a non- 
linear relativistic Fokker-Planck dynamics which incorporates non-extensive 
statistics and anomalous diffusion. 

A class of anomalous diffusions are currently described through the non- 
linear Fokker-Planck equation (NLFPE) 



J(y)[f(y,t)r + D^[f(y,t)r 



(7) 



where D and_J are the diffusion and drift coefficients, respectively. Tsallis 
and Bukman Lj have shown that, for linear drift, the time dependent solution 
of the above equation is a Tsallis-like distribution with q = 1 + ji — v. The 
norm of the distribution is conserved at all times only if fi = 1, therefore we 
will limit the discussion to the case v = 2 — q. 

Imposing the validity of the Einstein relation for Brownian particles, we 
can generalize to the relativistic case the standard expressions of diffusion and 
drift coefficients as follows 

D = aT , J(y) = a m± sinh(y) = ap\\ , (8) 

where p\\ is the longitudinal momentum and a is a common constant. It 
is easy to see that the above coefficients give us the Boltzmann stationary 
distribution in the linear Fokker-Planck equation (q = v = 1) (such_a result 
cannot be obtained if one assumes a linear drift coefficient as in Ref.El) while 
the stationary solution of the NLFPE (Q) with v — 2 — q is a Tsallis-like 
distribution with the relativistic energy E — mj_ cosh(y): 

r 1/(1-9) 
f q (y,mi_) = jl - (1 -q)Pm±cosii(y)j . (9) 
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Basic assumption of our analysis is that the rapidity distribution is not 
appreciably influenced by transverse dynamics, which is considered in thermal 
equilibrium. This hypothesis is well confirmed by the experimental datamLj 
and adopted in many theoretical worksO. Therefore, the time dependent ra- 
pidity distribution can be obtained, first, by means of numerical integration 
of Eq. (0) with initial S- function condition depending on the value of the ex- 
perimental projectile rapidities and, second, by integrating such a result over 
the transverse mass m± (or transverse momentum) as follows 

dN f°° 

-j-(y,t) = c / cosh(y) f q (y,m±,t)dm± , (10) 

"•V Jm 

where c is the normalization constant fixed by the total number of the parti- 
cles. The calculated rapidity spectra will ultimately depend on the nonexten- 
sive parameter q only, since there exists only on£. "interaction time" r = at 
which reproduces the experimental distribution EJ. 




Figure 1. Rapidity spectra for net proton production (p — p) in central Pb+Pb collisions at 
158A GeV/c (grey circles are data reflected about y cm = 0). Full line corresponds to our 
results by using a non-linear evolution equation (q = 1.25), dashed line corresponds to the 
linear case (q = 1). 



In Fig.l we show the calculated rapidity spectra of net proton compared 

S'th the experimental NA49 data from central Pb+Pb collisions at 158 GeV/c 
The obtained spectra are normalized tp_164 protons and the beam rapidity 
is fixed to y cm = 2.9 (in the cm. frame) Ej. The full line corresponds to the 
NLFPE solution (|ic| ) at r = 0.82 and q — 1.25; the dashed line corresponds to 
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the solution of the linear case (q = 1) at r = 1.2. Only in the non-linear case 
(q ^ 1) there exists a (finite) time for which the obtained rapidity spectra 
well reproduces the broad experimental shape. A value of q =/= 1 implies 
anomalous supxxdiffusion in the rapidity space, i.e., [y(t) — yM(t)] 2 scales like 
t a with a > 1 III 



5 Conclusions 



The nonextensive statistics appears suitable to evaluate physical observables 
recently measured in heavy ion collision experiments. The physical motivation 
for the validity of a non-conventional statistical behavior can be related to the 
presence of memory effects and long range interactions at the early stage of 
the collisions, even if a microscopic justification of these effects is still lacking. 
A rigorous determination of the conditions that produce a nonextensive statis- 
tical regime should be based on microscopic calculations relative to the parton 
plasma originated during the high energy collisions. Non-perturbative QCD 
effects in the proximity of hadronic deconfincmcnt could play a crucial role in 
the determination of the quantum kinetic evolution of the system toward the 
equilibrium a. 
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